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I.  RESULTS  FROM  PRIOR  SUPPORT 


1.  Simulation  Study  of  a  New  Mechanism  for  Excitation  of  Kinetic  Waves  in  a 
Magnetoplasma 

Nishikawa  et  al.  [1988]  have  investigated  the  new  ion-cyclotron-like  waves  by  a  localized 
transverse  electric  field  by  means  of  simulation  with  the  assistance  of  the  nonlocal  kinetic 
theory.  The  linear  theory  shows  that  the  growth  rates  of  the  kinetic  Kelvin- Helmholtz  (K-H) 
modes  are  strongly  reduced  with  increasing  u  =  k\\/ky,  and  they  become  unstable  only  when 
6  =  kyp?  <  1  and  kyL  ss  1,  where  L  is  the  scale  length  associated  with  the  transverse  electric 
field.  On  the  other  hand,  the  new  modes  have  larger  frequencies  and  become  unstable  at  larger 
b  >  1  and  kyL  >  1  [Ganguli  et  al.,  1988]. 


Simulation  results  show  that  ion-cyclotron-like  waves  are  excited  in  regions  where  E  X  B 
drift  is  localized.  Linear  growth  rates  of  several  modes  are  obtained  from  the  wave  analysis 
of  the  simulation.  This  linear  analysis  shows  that  the  (0,  4)  mode  corresponds  to  large  6,  and 
large  kyL  has  the  maximum  growth  rate.  Clearly,  these  are  not  K-H  modes.  Further,  the 
simulation  results  show  that  density  gradients  help  to  enhance  tne  growth  rates.  However,  like 


the  K-H  mode,  the  real  frequencies  of  this  instability  are  approximately  proportional  to  kyV g, 
where  kg  is  the  peak  value  of  the  E  x  B  drift. 

Nonlinear  phenomena  such  as  diffusion  and  coalescence  of  vortices  are  investigated.  In  the 
linear  stage  smaller  vortices  are  generated  and  larger  vortices  with  the  lower  frequencies  are 
dominant  in  the  nonlinear  stage.  In  the  nonlinear  stage  ions  diffuse  strongly  due  to  large-scale 


□  0 


Recently,  we  have  investigated  the  electrostatic  waves  driven  by  the  combined  effects  of  a 
localized  transverse  electric  field  and  parallel  electron  drifts  by  means  of  simulation  with  the 
assistance  of  the  nonlocal  kinetic  theory  [Nishikawa  et  al.,  1989b, c,d]. 

We  have  performed  a  number  of  simulations  for  this  instability.  Simulation  results  show 
that  electrostatic  waves  are  excited  in  the  regions  where  the  ExB  drift  is  localized  in  the 
simultaneous  presence  of  parallel  electron  drifts  and  transverse  electric  fields.  Simulations  with 
only  parallel  electron  drifts  or  transverse  electric  fields  show  that  no  instability  grows  out  of 
the  thermal  noise.  The  simulation  with  both  the  parallel  electron  drift  and  the  transverse 
electric  field  shows  the  growing  waves  out  of  the  thermal  noise.  The  Doppler  shift  due  to  the 
ExB  drift  can  lower  the  phase  velocity  of  waves  along  the  magnetic  field.  Then  this  makes 
wave-particle  resonance  possible  for  smaller  V£,  which  leads  to  this  instability. 

The  nonlinear  phenomena  such  as  diffusion  and  coalescence  of  vortices  are  investigated.  In 
the  linear  stage,  smaller  vortices  are  generated  and  larger  vortices  with  the  lower  frequencies 
are  dominant  in  the  nonlinear  stage.  In  the  nonlinear  stage  the  ions  diffuse  strongly  due  to 
the  large  scale  vortices. 
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2.  Beam  Instability  in  the  Foreshock 


As  an  application  of  the  simulation  method  used  in  the  proposed  research  (Broadband 
electrostatic  noise),  the  beam  instability  in  the  foreshock  has  been  investigated.  Electrons 
backstreaming  into  the  Earth’s  foreshock  generate  waves  near  the  plasma  frequency  by  the 
beam  instability.  Two  versions  of  the  beam  instability  exist:  the  ‘reactive’  version  in  which 
narrowband  waves  grow  by  bunching  the  electrons  in  space,  and  the  ‘kinetic’  version  in  which 
broadband  growth  occurs  by  a  maser  mechanism  [Cairns,  1887a,  b,  and  references  therein].  Re¬ 
cently,  Cairns  [1987b]  has  suggested  that  (1)  the  backstreaming  electrons  have  steep-sided  ‘cut¬ 
off’  distributions  which  are  initially  unstable  to  the  reactive  instability,  (2)  the  back-reaction 
to  the  wave  growth  causes  the  instability  to  pass  into  its  kinetic  phase,  and  (3)  the  kinetic 
instability  saturates  by  quasilinear  relaxation. 

Cairns  and  Nishikawa  [1989]  have  performed  two-dimensional  simulations  of  the  reactive 
instability  for  Maxwellian  beams  and  cutoff  distributions.  The  results  of  the  simulations  are 
consistent  with  suggestions  (1)  and  (2)  above.  In  addition,  we  have  demonstrated  that  the 
reactive  instability  is  a  bunching  instability,  and  the  reactive  instability  saturates  and  passes 
over  into  the  kinetic  phase  by  particle  trapping.  We  found  that  the  kinetic  growth  occurring 
after  saturation  of  the  reactive  instability  is  presumably  due  to  the  spatially  localized  gradients 
in  y—  V||  phase  space.  Both  simulation  results  and  numerical  solutions  of  the  dispersion  equation 
indicate  that  the  center  frequency  of  the  intense  narrowband  waves  near  the  foreshock  boundary 
may  be  between  0.9wpe  and  0.98u>pe,  rather  than  being  above  upe  as  previously  believed. 
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3.  Whistler  Mode  Driven  by  the  Spacelab  2  Electron  Beam 


During  the  Spacelab  2  mission  while  an  electron  beam  was  being  ejected  from  the  shuttle, 
the  Plasma  Diagnostics  Package  (PDP)  detected  a  clear  funnel-shaped  emission  that  is  believed 
to  be  caused  by  whistler-mode  emission  from  the  electron  beam  [Gurnett  et  al.,  1986].  In  order 
to  understand  the  mechanism  of  this  emission,  the  simulations  have  been  performed  using 
a  three-dimensional  magnetostatic  code  for  low-/?  plasmas  in  which  the  beam  electrons  are 
initially  located  in  the  column  [Nishikawa  et  al.,  1989a].  In  order  to  simulate  the  continuous 
electron  ejection  from  the  shuttle,  the  simulations  were  also  performed  with  the  recycling  of 
the  beam  electrons.  The  beam  electrons  excite  whistler  waves  and  lower  hybrid  waves.  The 
brief  fluid  theory  based  on  the  magnetostatic  code  was  checked  with  tne  simulation  results. 
The  propagating  whistler  mode  was  identified  with  the  theory.  The  simulation  results  show 
that  the  quasi  perpendicularly  (the  angle  between  the  magnetic  field  and  the  wavenumber  is 
larger  than  50°)  propagating  whistler  waves  have  larger  amplitude  whose  real  frequencies  are 
smaller  than  the  local  electron  cyclotron  frequency.  This  fact  is  consistent  with  the  fact  that 
the  funnel-shaped  emission  is  observed  below  the  electron  cyclotron  frequency  away  from  the 
beam  electron.  The  beam  electrons  initially  in  the  column  diffuse  radially  as  well  as  slow  down 
due  to  the  E  x  B  caused  by  the  excited  waves.  The  acceleration  of  the  beam  electrons  also 
takes  place  due  to  the  excited  whistler  waves. 

In  order  to  compare  with  the  PDP  data,  the  local  magnetic  fields  Z?xy  and  the  perturbed 
electric  fields  EVtZ  are  diagnosed  at  the  several  points  in  the  simulation  system.  The  results 
show  that  the  waves  are  radially  excited  by  the  beam  electrons  localized  in  the  center  of  the 
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system.  The  wave  spectra  of  the  electric  fields  Ey>z  diagnosed  at  x  =  31Ae,  y  =  16Ae,  and 
z  =  17Ae  show  the  several  kinds  of  waves  generated  by  the  beam  electrons.  The  analysis 
shows  that  the  lower  hybrid  waves  and  whistler  waves  are  excited.  The  frequency  range  of 
these  spectra  extends  from  up{  to  beyond  upe  which  is  in  qualitative  agreement  with  the 
PDP  data.  As  the  PDP  data  show,  the  intense  electrostatic  narrowband  emission  around  the 
electron  plasma  frequency  has  been  observed  dominantly  in  the  spectra  of  the  electric  field  ( Ez ) 
with  and  without  the  recycling  of  the  beam  electrons.  This  wave  is  identified  as  the  parallel 
(quasi  parallel)  whistler  wave.  This  simulation  result  is  in  good  agreement  with  the  PDP  data 
(Gurnett  et  al.,  1986].  The  spectra  of  the  magnetic  fields  diagnosed  at  the  same  position  show 
that  whistler  modes  are  excited.  The  cB/ E  was  calculated  for  the  whistler  mode.  In  the  case 
without  the  recycling  of  the  beam  electrons,  the  cB/E  is  approximately  0.35  for  the  whistler 
mode  (around  0.368up<.),  which  is  smaller  than  that  calculated  from  the  PDP  data  [Gurnett 
et  al.,  1986]. 
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4.  Simulation  of  Electron  Cyclotron  Harmonic  Waves  by  AMPTE  Lithium  Release 


We  have  studied  the  generating  electrostatic  turbulence  in  the  lithium  cloud  due  to  solar 
wind  particles  which  are  flowing  through  it  [Nishikawa  et  al. ,  1989e].  Several  simulations  have 
been  performed  in  order  to  understand  the  processes  generating  electron  cyclotron  harmonic 
waves. 

The  simulations  show  that  electron  cyclotron  harmonic  waves  are  excited  by  mainly  the 
solar  wind  electrons.  The  spectra  of  perturbations  of  electrostatic  potential  exhibit  electron 
cyclotron  harmonic  waves. 

The  overall  frequency  range  is  obtained  by  observing  the  electric  fields  at  some  points  in 
the  simulation  region.  Power  spectra  of  the  electric  fields  exhibit  electron  cyclotron  frequency 
and  higher  harmonics  which  is  simular  to  the  spectrum  obtained  by  the  satellite. 

Possible  mechanisms  of  generating  electron  cyclotron  harmonic  waves  have  been  discussed 
in  the  previous  work  [Roeder  et  al.,  1987].  In  this  simulation  study,  the  possible  energy 
sources  of  the  instability  are  included  automatically.  The  relative  velocities  among  different 
species  become  the  energy  source  of  instabilities.  The  relative  velocity  between  solar  wind 
electrons  and  lithium  ions  would  be  one  possible  energy  source  for  the  excitation  of  ECH 
waves.  The  fact  that  the  instability  can  be  excited  without  new-born  cold  electrons  support 
this  mechanism.  Theoretically,  the  mode  couplings  among  harmonics  occur  mainly  due  to  the 
electron  perpendicular  drift,  which  lead  to  the  instability  [Tataronis  and  Crawford,  1970]. 

According  to  Figure  1  by  Gurnett  et  al.  [1986b],  the  higher  frequency  ranges  of  waves  look 
very  simular  at  both  cases  with  ntl  ss  npr  and  rq,  <C  npT.  In  this  study,  we  found  that  ECH 
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waves  are  excited  also  in  the  case  of  n;,-  =  npr.  The  satellite  data  show  that  ECH  waves  can  be 
observed  only  in  some  range  of  lithium  ion  density  (0.03  -  0.18  cm-3)  associated  with  lithium 
drift  energy  [Roeder  et  al.,  1987].  The  observations  show  that  in  the  case  of  nu  as  npr  ECH 
waves  are  not  found.  The  simulation  results  show  that  in  the  same  case  of  nu  —  npr  ECH  waves 
are  excited.  This  discrepancy  may  be  explained  by  the  setting  of  simulation  geometry  due  to 
the  two-dimensional  code.  If  the  simulation  plane  is  set  nearly  parallel  to  the  magnetic  fields 
(B0y/ Bo  <  1),  the  angle  between  the  magnetic  fields  and  the  drift  velocity  of  hot  electrons  and 
protons  becomes  very  small.  Then,  in  the  case  of  nu  =  npr,  the  parallel  drift  of  solar  wind 
protons  and  electrons  may  excite  ion  acoustic  waves  more  dominantly  than  ECH  waves,  which 
is  not  included  in  this  study  [Gurnett  et  al .,  1986b;  Ma  et  al.,  1987]. 

Further  study  is  required  on  the  observed  plasma  velocity  distributions  and  comparisons 
of  simulation  wave  results  with  the  wave  power  spectra.  More  simulations  are  required  for 
further  progress  in  the  interpretation  of  electrostatic  noise. 
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